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Abstract 

We study supersymmetric models with a GUT-sized extra dimension, where both 
the Higgs fields and the SUSY breaking hidden sector are localized on a 4D brane. 
Exponential wave function profiles of the matter fields give rise to hierarchical struc- 
tures in the Yukawa couplings and soft terms. Such structures can naturally explain 
hierarchical fermion masses and mixings, while at the same time alleviating the 
supersymmetric flavour problem. We discuss two sources of supersymmetry break- 
ing, radion mediation and brane fields, and perform a detailed numerical analysis, 
thoroughly taking into account the proliferation of unknown 0(1) coefficients that 
occurs in this class of models. It turns out that additional assumptions on supersym- 
metry breaking are necessary to evade the stringent experimental bounds on lepton 
flavour violation. The favourable regions of parameter space arc then examined with 
regards to their LHC phenomenology. They generically feature heavy gluinos and 
squarks beyond current bounds. Lepton flavour violation in SUSY cascade decays 
can give interesting signatures. 



1 Introduction 



In the Standard Model, the three generations of quarks and leptons follow a peculiar 
pattern of hierarchical masses and mixings. Localizing the Standard Model matter fields 
in the bulk of a compact extra dimension, for instance on a slice of AdSs naturally 
leads to such flavour hierarchies (2||3]. 

For a strongly warped extra dimension, warping could explain the discrepancy be- 
tween the Planck scale and the electroweak scale fl]. In that case the lowest massive 
Kaluza-Klein modes should have masses around a TeV. In less strongly warped (or even 
unwarped) models, the KK scale can be only a few orders of magnitude below the Planck 
scale, and the electroweak hierarchy problem can be solved by TeV-scale supersymme- 
try. Wave function localization still accounts for the Yukawa hierarchies, and the warped 
internal space still allows for an interpretation of the 5D model as a holographic dual of 
some strongly coupled CFT [4j. If the KK scale is not too far from the 4D GUT scale, 
such a model may even be compatible with standard (logarithmic) MSSM gauge cou- 
pling unification, with the GUT group broken at the compactification scale by boundary 
conditions. 

As an additional benefit of such models, localizing matter fields in 5D can allevi- 
ate the SUSY flavour problem (see, for example, jsjjsj). This is what happens if the 
Higgs fields are brane fields, and supersymmetry breaking is localized on the same brane 
(and possibly in the gravitational background). In that case the trilinear soft terms 
follow a hierarchy structure similar to that of the the Yukawa couplings, and will thus 
be approximately diagonal in the fermion mass eigenstate basis. Off-diagonal scalar soft 
masses may also be suppressed. A related mechanism was originally advocated in a 4D 
setup, with the visible sector fields acquiring large anomalous dimensions from their cou- 
plings to a strongly coupled CFT [Oj . These two pictures can be argued to be related by 
AdS / CFT duality [5j . Whether or not FCNCs in such a setup are sufficiently suppressed 
to evade experimental bounds is however heavily model-dependent. 

This question is the subject of the present paper. We will study, using a concrete 
example model as our benchmark, to what extent wave-function localization in 5D is 
enough to suppress flavour- violating processes, and to what extent some additional mech- 
anism or some residual tuning is still needed. We will discuss mass spectra and collider 
signatures in the favoured regions of parameter space. In the present work we focus on 
FCNCs and neglect CP violation. 

Previous studies of this subject have revealed that wave-function suppression alone 
is in general not enough to evade the stringent experimental bounds, and that some 
additional mechanism or residual tuning of the parameters is needed. The authors 
of [5] suggested that 5D mass terms, which determine the localization properties of the 
bulk zero modes, might be required to be quantized in discrete units by the underlying 
fundamental theory. With suitably chosen discrete 5D masses, and some moderate 
residual tuning, they argued that all flavour constraints could be satisfied. In [T] a 
U(l) symmetry was used to forbid some of the couplings between the SUSY breaking 
sector and the visible sector, thus enforcing flavour- diagonal a-terms. Flavour violation 
in the scalar soft masses was then argued to typically be sufficiently suppressed by wave- 
function localization alone. A recent study of flavour violation in a 5D warped model 
with gaugino mediation was conducted in [lOj. 

Much of our analysis is sufficiently broad to be representative for any 5D SUSY GUT 
with the Higgs and SUSY breaking sectors localized on the same brane. In particular 
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the results on the low-energy spectrum should generalize, at least qualitatively, to any 
such model (for instance, those constructed e.g. in [7|). However, when we do need to 
work with a concrete model for definiteness, we choose the "holographic GUT" model 
of Nomura, Poland and Tweedie (NPT) [6]. In the NPT model, there is a warped extra 
dimension, and the bulk gauge group is SU(6). It is broken by boundary conditions 
to SU(5) X U(l) on the UV brane, and by the VEV of an adjoint brane field S to 
SU(4) X SU(2) X U(l) on the IR brane. This gives essentially the Standard Model 
gauge group in the 4D effective field theory. Matter fields are localized in the bulk, and 
boundary conditions are chosen such that their zero modes furnish precisely the matter 
content of the MSSM. The MSSM Higgs fields are pseudo-Goldstone bosons arising 
from S [11,12 , and the symmetry breaking structure results in a GUT-scale degenerate 



Higgs mass matrix 13 , thus automatically solving the fi problem. Supersymmetry 
is broken both by some hidden sector chiral superfields on the IR brane, and by the 
5D gravitational background (i.e. by non-vanishing F-terms of the radion and chiral 
compensator superfields in 4D language). Models of this type have been shown to be able 
to give realistic low-energy mass spectra when assuming flavour-blind supersymmetry 
breaking [13| . Therefore it is interesting to now relax the assumption of flavour-blindness 
and study the impact on flavour observables. 

An important aspect of our work is that we carry out an extensive numerical analysis, 
giving quantitative results for flavour constraints and collider phenomenology. In that 
respect we go beyond the previous literature, in particular jsjjs], which contain only 
qualitative discussions. 

To avoid confusion it may be worth noting that our study is conceptually unrelated 
to work on flavour violation in Randall-Sundrum modelsj^ Such models are typically 
non-supersymmetric, since the hierarchy problem is solved by the strongly warped fifth 
dimension. However, with the Standard Model fermions in the bulk, KK mode ex- 
change represents a potentially dangerous source for FCNCs. In the models which we 
are concerned with, by contrast, the KK modes are far too heavy to contribute to flavour 
violation. Flavour- violating processes can instead be mediated by the exchange of MSSM 
superpartners. 

This work is organized as follows: We start by reviewing first the supersymmetric 
flavour problem in Section [2| and then the generation of hierarchical fermion masses and 
mixings from wave-function localization in Section |3j We present our parameterization 
of supersymmetry breaking in warped 5D models in Section [4} Estimates for the mag- 



nitudes of soft parameters from naive dimensional analysis are presented in 4.1, for the 
case that SUSY breaking mediation is dominated by brane fields. Soft terms for the 
case that radion mediation dominates are given in |4.2[ In Section [5] we explain details 
about the numerical analysis which we perform to scan over a large number of models, 
with the results presented in Section [6} Conclusions are given in Section [Tj Appendix [A] 
contains a brief review of an example model in which SUSY is broken and the radius 
of the extra dimension is stabilized, both of which we assume to be the case in the 
main text without providing an actual mechanism. Appendix |B] contains some analytic 
expressions for the diagonalization of hierarchical Yukawa matrices, which were used in 
the numerical analysis. 



^See, for instance, |14j and references therein. 
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2 The SUSY flavour problem 



We start with a brief recapitulation of the flavour problem in supersymmetry, see e.g. 15 , 
to set up our notation and terminology. The Lagrangian for the flavour sector of the 
R-parity symmetric MSSM reads 

+ J cfe {y^j HuQiUj + yf^ H^QrOj + yf^ HdUE,) + h.c. 

+ ml., qjqj + ufn^ + dldj + mf^jljlj + ml^jelej 
+ {afj huqiUj + afj hdqidj + afj hdhej + h.c.) . 

Here Q, U and D are the quark superfields and E and L are the lepton superfields (we 
are omitting right-handed neutrinos) with q, u, d, I and e their scalar components. 
and Hd are the Higgs superfields and and hd their scalar components. The couplings 
y^ , y", y^ , m^, m^, m^, ml, rn^, , a° and are complex 3x3 matrices in flavour 
space; the soft mass matrices with X = U, D, Q, L, E are restricted to be hermitian. 

Flavour rotations X — )■ UxX, with X = any matter superfield and Ux a unitary 
3x3 matrix, leave the kinetic terms invariant but change the couplings. By appropriately 
rotating E — t- UeE and L — )■ UlL, the Yukawa matrix y^ can be diagonalized. 

Likewise, it is possible to choose either or diagonal. However, either of these 
two choices fixes Uq up to a phase, so and y^ cannot be chosen diagonal simultane- 
ously: There is CKM mixing in the quark sector of the Standard Model. It is of course 
possible to go to a field basis where S\J{2)l is non-linearly realized, by splitting the weak 
doublet according to Q = {Ul, Dl)^ (in a gauge where the Higgs expectation values will 
eventually be (hy) = (0, v")^ and (hd) = {vd, 0)^). One can then perform independent 
flavour rotations on Ul and D l and thus diagonalize both up- type and down- type mass 
matrices simultaneously; this is the "super-CKM basis", in which the quarks are mass 
eigenstates but no longer weak interaction eigenstates. 

In the limit of vanishing soft terms, flavour-changing processes are suppressed in the 
MSSM as they are in the Standard Model. In particular, there are no FCNCs at the 
tree level. Processes such as neutral meson mixing K ■v^ K and D ^ D are present at 
one loop, but by the GIM mechanism the corresponding box diagrams are suppressed 
by {ml — m\)/m^r or {ml — m'^)/m^^ (where mu,d,c,s are the quark masses, which of 
course coincide with the squark masses in the supersymmetric limit). 

By contrast, in general the soft terms and will be non-diagonal in the CKM 
basis. Furthermore, in general the squark masses for the first two generations have no 
reason to be small or near-degenerate, so the GIM suppression is lost. Generic soft 
mass matrices and trilinear terms will therefore give unacceptably large contributions 
to strongly constrained processes such as flavour- violating lepton or meson decays, or 
neutral meson mixing. 

This problem should be addressed by imposing some specific structure on the soft 
masses and trilinear soft termsj^ For instance, if the a-terms a^, a^, are proportional 

^Of course the sfermions might as well be very heavy, some tens of TeV at least, such that they 
effectively decouple. This scenario, however, is disfavoured by naturalness arguments, since for SUSY to 
provide a natural solution to the electroweak hierarchy problem, the superpartner masses should not be 
too far above the electroweak scale. 
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to the respective Yukawa matrices, then they will evidently be diagonal in the CKM 
basis. Furthermore, if ~ 1 for all X, then the soft mass matrices are unaffected 
by flavour rotations and remain diagonal in the CKM basis, with degenerate entries. 
Exact or approximate patterns like these can be the consequence of some flavour-blind 
mechanism of supersymmetry breaking mediation, or they can result from horizontal 
symmetries, or from wave-function localization in an extra dimension. In this work we 
are investigating the latter mechanism. 



3 Flavour hierarchies from locahzation 

We will now review how localization of matter fields in an extra dimension can naturally 
generate large hierarchies in fermion masses and mixings, and how it may at the same 
time ameliorate the supersymmetric flavour problem. 

To be precise, we will consider a 5D model, compactified on an interval such that 
the KK mass scale is close to the 4D GUT scale Mqut ~ 2 • 10^^ GeV. The zero modes 
contain the MSSM fields, with chiral matter obtained from 5D bulk hypermultiplets. The 
5D SUSY Lagrangian allows for hypermultiplet mass terms, which will not decouple the 
zero modes, but rather distort their wave function profiles in the fifth dimension. The 
resulting zero mode wave functions are exponentially localized towards one of the branes. 

The simplest example is given by a fiat extra dimension obtained as the Z2 orbifold 
of a circle of radius 27ri2. A 5D hypermultiplet can be decomposed into two 4D chiral 



superfields {X, X'^), whose action is 16 



S = j (fx j^^ dy j d^9 (^X^X + X^X^^^ + j d^e {X^dyX + MXX^)+ h.c. 



If X'^ is odd and X is even under the orbifold projection, X will have a zero mode whose 
profile is given by 

xo{y) ~ e-^'y. (3) 

Similarly, in a warped extra dimension with metric ds'^ = e~^ydx'^ + dy"^, the action 
for a hypermultiplet reads flTj 



S = j d^x £^ dy j d^9 e-^^y [x^ X + X^X^^ 

+ j d^9e-^^y (^X^dyX - ^XdyX^ + MXX''^ + h.c 

Defining c = M/k, as is common convention, the X scalar zero mode has a profile [3] 

xo{y) = -^e^l-^>^y. (5) 



(4) 



For later convenience, we also define the profile factor / to be xo{y) evaluated on the 
y = ttR brane and with one power of the warp factor absorbed, 



/-kR 



(6) 
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After integrating over y one obtains for the zero mode kinetic action in four dimensions 



S^ = J (fx J (fOYXlXQ, (7) 

where 

j\-c)2-KkR _ 1 

Y=vi ^ -■ (8) 

[\-c) 2TrkR ^ ' 

For c > I the zero mode is locahzed towards the y = ttR (IR) brane, whereas for c < | 
it is locahzed towards the y = (UV) brane. For c = ^ it has a flat profile. In the 
following we always assume c > — ^, because for c < — ^ some massive modes with 
twisted boundary conditions could become exponentially light. With c > —5, all KK 
excitations decouple around the scale ke~'"^^, which we choose around Mqut- 

One may also introduce purely 4D chiral superfields which are entirely localized on 
one of the branes. For an IR-brane Higgs superfield H the action can contain, besides 
the kinetic terms, also Yukawa couplings between H and the even 4D chiral parts Xi of 
5D hypermultiplets {Xi,Xf): 

5brane = j (fx j dy j (fOe-^^yH^H + j (fO e'^^y hij HX.Xj + h.c. 6{y-irR) (9) 

The effective 4D Yukawa coupling for canonically normalized zero modes reads 

Uii = hij "^^ "^'^ ,„ (no summation) , (10) 

(YiY,y/2 

where the Yj are defined analogously as in Eq. ([s]), and the fi as in Eq. ([6]). Since the 
fi depend exponentially on the Cj, hierarchical Yukawa matrices are obtained from 0(1) 
c-parameters and anarchical hij coefficients. 

As a concrete example, let us review the flavour sector of the holographic GUT 
model of [6]. The gauge symmetry in the bulk is SU(6), broken to SU(5) x U(l) by 
boundary conditions on the UV brane. The bulk zero modes are contained in three 4D 
chiral superflelds Ti (each containing a 5 of SU(5)), three chiral superfields % (each 
containing a 10), and three chiral superfields J\fi (each containing a singlet). These 
zero modes are identified with the MSSM matter fields. On the IR brane the boundary 
conditions preserve the SU(6) bulk gauge symmetry. However, there is an IR brane 
superfield S in the adjoint of SU(6) whose expectation value breaks SU(6) spontaneously 
to SU(4) X SU(2) X U(l). The low-energy gauge group is given by the intersection of 
SU(5) X U(l) and SU(4) x SU(2) x U(l) in SU(6), which is the Standard Model gauge 
group (apart from an extra U(l) which is Higgsed on the UV brane). Two weak doublet 
components of S remain massless, and are identified with the MSSM Higgs doublets. 
The other S components either acquire supersymmetric masses directly from the scalar 
potential, or are eaten by the Higgs mechanism. The field content is listed in Table [l} 

A possible choice for the scales in this model is to fix the cutoff scale of the 5D theory 
to be around the 4D reduced Planck scale, = 2 • 10^® GeV say, and the KK scale 
k e-""^^ to be slightly below the GUT scale, k e-""^^ = 10^^ GeV. The AdS curvature k 
is chosen to lie in between, k = 2 • 10^^ GeV, and kR = 1. This gives a somewhat large 
size of the extra dimension in units of the cutoff, around ttRM^, = 30. 

The Yukawa couplings arise from 

CD6{y- ttR) j (fe [hjj i:TiTj + hfj T.TiFj + h^j ^FiMj) + h.c. (11) 
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field 


SU(6) rep. 


massless mode 


localization 


TTC 


70 
70 


5 


bulk 
bulk 


Ti 


20 
20 


10 


bulk 
bulk 


Mi 


56 
56 


1 


bulk 
bulk 




35 


2l/2 © 2_i/2 

of SU(2) X U(l) 


IR brane 



Table 1: Field content of the NPT model, listing the SU(6) representations and the massless 
modes' SU(5) representations (if applicable). 



The IR brane couplings /ijj can be estimated using naive dimensional analysis; their 
typical magnitude up to 0(\) uncertainty is 67r^/M*. 

For illustration consider the following choice of c-parameters: Setting 

CTi ~ 2' ^T ~ « cj-3 « -, CTg « CA^i ^ cm^ ~ CA/-3 - (12) 

we obtain in the quark sector 

V" = >^ifl A|e , y°=\ A^e^ A|62 , (13) 

\ ^^13^ -^23^ ^^33 / \ -^31^ ^^32^ -^33^ / 

where e w e'""^^ w ^, and where the X"^-" are of the order X'^f = O (Qir/M^R) ^ 0(1). 
Note that is symmetric, because it arises from the symmetric 'ETiTj coupling in 



Eq. (11). In the lepton sector one has 



^11^ ^12^ -^13 \ / Afl '^f2 ^13 

= I A|62 Xl,e Af3 , = A| A^^ X^, | . (14) 

^31*= AfgE Af3 / \ X^i X^2 A33 

Hierarchical Yukawa matrices of this type are well-known to roughly give the observed 
masses and mixings. Unwanted SU(5) relations can be avoided by taking into account 
contributions from higher-dimensional operators, in particular higher powers of S. 



The assignment Eq. (12) can be refined to even better reproduce the known fermion 
masses and CKM angles. We will describe in detail how we fit the c-parameters and the 
Xij in Section [5| 



4 Supersymmetry breaking 

SUSY breaking can be parameterized by F- and D-type spurious in the Kiihler potential, 
which we generally denote by $ and A, and by F-type spurious in the superpotential 
denoted by A. Omitting the neutrinos from now on, the 4D Lagrangian can be brought 
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into the form 



C= d 



Y//„ + i^Hj'' + h.c. ) + ^ny] HiHu + (Hu o Ha) 



+ / d'e 



+ ( Y,^; + ( + h.c.) + uJUj + ([/ o {D, Q, L, E}) 

u q2\ tt tt I ( ;r.D , 7d q2 



(15) 



+ h.c. 



Here the Y, A, $ and A are c- number functions of the compactification radius and 
of the expectation values of hidden sector fields. In particular, in the absence of brane 
kinetic terms (which are generically subdominant at large volume), the wave-function 
coefficient Y for the bulk fields is as in Eq. (Isl): 



Y^. - 6- ■ - - 



,2 cf)2TTkR ' 

where X = U, D, Q, L, E. 

Again for X = any matter field, we define rescaled quantities by 







y/2 




a jj 













-'-^ 1/, 



Y'^'^Y'^.Yf/ 



1/2'' 



A^. 



13 



A^- 



ii jj 



1/2' 



A^. 



The Yukawa matrices for canonically normalized fields are then 



1/2 • 



1/2 ■ 



Y^Y'^Yff 

II 33 



\ II 33 

The scalar soft masses for matter fields are 



1/2^ 



(y|IY|^.Yh,) 



1/2 



($^)t ^x _ 



and the trilinear terms are given by 



a" 


= ($« 


)V + y" 


(*T + */^.y" 


- A^ 


a" 


= (*^^ 


) V + 


(*-)^ + *^,,y- 


- A^ 




= m 


V + y"( 







(16) 



(17) 



(18) 



(19) 



(20) 



There are two natural possibilities for the origin of supersymmetry breaking in our 
5D setup: Supersymmetry breaking by the radion multiplet [18] , or supersymmetry 
breaking by additional brane fields. In radion-mediated SUSY breaking, the radion 
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superfield (whose lowest component sets the radius of the extra dimension) acquires 
an F-teim expectation value. This scenario is quite predictive since all the couplings 
are essentially determined from geometry. By contrast, if there are additional SUSY- 
breaking fields Zj on the branes, their couplings to the visible sector are additional 
free parameters. The general case will be a mixture of these two; see the Appendix 
for a concrete model in which both brane field SUSY breaking and radion mediation 
contribute to the soft terms. 

In order to alleviate the flavour problem, we make the crucial assumptions that the 
MSSM Higgs fields are brane fields (as is the case e.g. in the NPT model) and that the 
Zj are localized on the same brane. With the profile functions of Eq. ([s]), and with the 
dependence on the radion multiplet restored [Tt], the 4D Lagrangian is 

g{|-cf )fc7r(T+T) _ ^ 

-^-c^)2TTkR 
+ {U^{D, Q, L, E}) 

+ I d^e (^(^e-'=-(^+^) (l + [^HAZi)e' + h.c.) + AhAZi)0') HiHu 

+ {Hu ^ Hd) (21) 

+ / d'^e ^^Yl [mZi)0' + h-c. + A^^{Zi)9'') uJUj 

ij 

+ {U^{D, Q, L, E}) 

+ j d^e VFbrane + h-C 

where 



J a-cr)27rkR 



ttR 

TkT 



-(c^+c^)7rfcT 

+ (/.,- + A'!^{Zj)e') ' ^ HdE.L, ' -3.fcT 



Here T = R+iB^+F'^ 9"^ +{feTm.ions) is the radion multiplet, with the fifth component 
of the graviphoton. ip = l+F^O"^ is the chiral compensator, an auxiliary chiral superfield 
whose F-term component F"^ is the scalar auxiliary of the 4D gravitational multiplet. 
SUSY breaking in the combined gravitational sector is parameterized by F'^ ^ and 
Ff 7^ 0. Note the non-standard dimensions of these fields: In our conventions T and ip 
have mass dimensions —1 and respectively, so F'^ and F'^ have mass dimension and 1. 
The spurious A and A capture the effects of SUSY breaking by the IR brane fields 
Zi (which we also take as a background). We have assumed negligible brane-kinetic 
terms for the bulk fields. 



The matter wave function normalization coefficients are as in Eq. (16), while the 
Higgs wave function normalization reads 



Yh„,, = e-^-'^^ . (23) 
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For the spurionic coefficients one obtains, using profile factors as defined in Eq. 



F'P - TTkF^ + $ 

g(|-cf)27rfciJ 



^ {\- cf)2TTkR 



-2-KkR 



1 ^e(|-cf)2-fci? ^ 



2R 



-I] J I 



e^2 



+ 



2 cf)27rkR 
1 



F'^ + fcvr 1 - - c- 



1 



(24) 



cf - c-)fcvr + + h.c. + A- ) /f // 



and 



3F^ - (cf + cf)k7TF'') hfj + K% 



-knR fD fQ 



3F^ - (cf + c^)A;7rF^) hf^ + Af, ) e''^^^ /f /, 



-kwR fE rL 
'j ■ 



(25) 



Here we have not hsted the spurions A//^^, which do not contribute to sfermion soft 
masses. The non-canonical Yukawa couplings ijij can be read off to be 





= e- 


-irkR jU jQ 


y?, 




-irkR jD jQ 






-TvkR jE jL 



(26) 



Using Eq. (18) and the wave function normalizations from Eqs. (16) and (23), this 



gives canonical Yukawa couplings as in Eq. (10) (note that the Higgs wave function 
normalization cancels the explicit warp factor). 

In order to discuss MSSM spectra, we also need gaugino masses. The 4D gauge field 
Lagrangian is 



4 ^ 

a 



1 ttT - 

+ ^ + nyzi)e' ] w^^w^ + h.c. + . . . 



(27) 



We have omitted terms irrelevant for gaugino masses, and neglected possible effects 



from bulk Chern-Simons terms 22,24 . In this expression a = 1, 2, 3 labels the Standard 



Model gauge factors, is the bulk gauge coupling, and the l/^uv t^rm originates from 
a gauge kinetic term on the UV brane. Note that this term is universal with respect 
to the Standard Model gauge fields, at least in models where the UV brane preserves 
SU(5). IR brane terms are irrelevant to our discussion and will be omitted. The gaugino 
masses are then given by 

1 o f-. ttF^' 



Ma 



-.94 



2 ' n» + 



9l 



(28) 



with l/g1 = l/^uv '^^/d'i- numerical analysis we take the gaugino masses to 

be equal, 

Mi = M2 = M3 = Ml/2 ■ (29) 

Universal gaugino masses are induced by the leading contributions in the NPT model, 
since in this model the GUT group is broken only by expectation values on the IR 
brane. Higher-dimensional operators involving powers of the GUT Higgs T, can lead to 
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gaugino mass splittings, but we assume that they are sufficiently suppressed. In a more 
general setup where the gauge symmetry is broken by boundary conditions on the SUSY 
breaking brane, the universality condition could be relaxed. 

It is evident that, in general, the SUSY breaking soft terms will not be flavour- 
preserving. If SUSY breaking is dominated by radion mediation, the scalar soft masses 
will be diagonal in the basis we have been using, but their eigenvalues will be non- 
degenerate; rotating to the CKM basis will therefore induce off-diagonal terms. The 
trilinear couplings, likewise, are not proportional to the Yukawa couplings. If SUSY 
breaking soft terms are induced predominantly by brane fields, there are even fewer 
constraints, since the brane spurious A, A are generally anarchic. 

However, a key assumption for this class of models was that the hierarchies in 
the fermion masses and mixings are originating mainly from wave-function localization. 
The same localization effects also leave their imprints the soft terms, and will induce 
similar hierarchies in the mass matrices and trilinears; such scenarios have been dubbed 



"flavourful supersymmetry" 19 . It is then reasonable to argue that the flavour problem 
should at least be alleviated, if not solved, in wave-function localization models. For 
instance, the basis changes used to switch to the CKM basis will approximately also 
diagonalize the trilinear terms, up to higher-order terms which are exponentially small. 
An extreme limiting case of this scenario would be to keep only the sfermion soft terms 
associated with the third generation, or even only those associated with the stop (and 
with additional flavour-blind contributions to the soft terms induced by RG running). 
The flavour constraints in this limiting case (which we previously studied in ^13j) are 
evidently far less severe. 

One of the aims of this work is to test the assertion that, moving away from this 
limiting case, wave-function localization still gives sizeable FCNC suppression, and thus 
substantially reduces the tuning required to get a realistic model. From the known 
Yukawa hierarchies we fix the 5D mass parameters cf . Using anarchical textures for the 
unknown couplings of the 5D theory, we then calculate the GUT-scale soft parameters, 
evolve them to the electroweak scale using their renormalization group equations, and 
calculate the resulting masses and mixings. These are finally compared to the existing 
bounds on fiavour observables. 



4.1 SUSY breaking dominated by brane fields 

In the case that the dominant source for the soft terms are brane fields Zj, there is a 
large number of unknown coupling parameters. If, in accordance with the holographic 
interpretation, the theory is taken to be strongly coupled on the IR brane and in the bulk, 
then the magnitude of these parameters can be estimated by naive dimensional analysis. 
We use loop factors £5 = 247r^ for 5D superfields (of dimension 3/2) and £4 = IGvr^ 
for 4D superfields (of dimension 1) [20]. Up to 0(1) uncertainty, the coefficient for a 
superpotential term composed of n bulk zero modes and m brane fields is 

(30) 



167r2 \^ j\/j3/2 J ^M,. 

while for a Kahler potential term it is 



(31) 
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Here M^, is again the cutoff scale, taken to be close to the reduced Planck scale in 4D. 
The 5D Yukawa couplings are thus estimated to be 



hi 



6^ 



(32) 



To estimate the IR brane spurion coefficients, we define a Goldstino superfield Z with 
expectation value (Z) = F^O"^ to be the combination of SUSY breaking brane fields 
which couples to the visible sector. Z will in general be composite, but we can treat 
it as a single, elementary, canonically normalized field for our purposes, absorbing any 
compositeness factors of Att or warp factors in the definition of . We obtain 



Gtt^F^ 



247r^ 



F' 



(33) 



247r3 F^ 
W' 



It is now convenient to define matrices 



"I] 



i j 



1/2 • 



(34) 



The nfj are dimensionless and hierarchical, with their largest entries at most of order 
of the top Yukawa coupling (as their structure is determined by the same profile factors 



as the Yukawa matrices). With Eqns. (34), (33), (24), (25), and (17) one obtains, up to 
0(1) uncertainty. 



F^ 



Hence the soft masses my 



A^ are dominated by A^ 



For the trilinear soft terms of Eq. ( 20 ) we find 



-F' 



where 



47ru^. 



(35) 



(36) 



(37) 



Clearly aV^ is dominated by the last two terms in Eq. (36). Provided that there are no 



accidental cancellations taking place between them, it is of the order 



(38) 



Analogous statements hold for the other a-terms. Comparing Eqs. (35) and (38), the 
largest a-terms will be around a factor of ~ larger than the largest scalar soft 



masses. 
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The NDA estimate for ff^ in Eq. ([27| is 



21 



67r2 



(39) 



with C a group-theoretical factor. For the NPT model with bulk gauge group SU(6) 
C2 (SU(6)) = 6. The canonically normalized gaugino masses of Eq. (28) are 



we set C 
thus 



Sn^gi F^ 



(40) 



V69I My 

The four-dimensional gauge coupling gets contributions from both the 5D bulk gauge 
coupling and a UV brane-localized kinetic term (cf. Eq. (|27|)): 



9l 



1 ttR 

5uv ai 



(41) 



The bulk gauge coupling can also be estimated from naive dimensional analysis; this 
gives 



9l 



C 

247r3 



(42) 



The UV brane, in the holographic picture, can be weakly coupled and thus g\jY can be 
smaller than its NDA value. In fact, unless R is rather large, (^uv must be small in order 
to obtain the proper unified gauge coupling g\ 0.5 from Eq. (41). 



Eqns. (42) and (40) lead to suppressed gaugino masses: 



Ma 



2_ 



0.05 



■El 



(43) 



Evidently, to have gaugino masses which are comparable with the other soft terms, 
IR brane terms cannot be the only source of SUSY breaking. Alternatively one could 
consider models which are weakly coupled also in the bulk and where, consequently, 
there is no large NDA suppression as in Eq. (42). 

Finally, we impose that the Higgs mass matrix be degenerate at the GUT scale: 



2 I I |2 



\B„ 



(44) 



This is the case in a large class of models 13 , in particular in the NPT model which we 
choose as our benchmark. The relation Eq. (44) has the advantage of constraining our 
free parameters somewhat more; it could however be relaxed in a more general setting. 
The NDA estimate for the Higgs mass parameters is 



\B„ 



IGvr^ 



F^ 



(45) 



4.2 SUSY breaking dominated by radion mediation 



It is instructive to also consider the opposite extreme case, where brane contributions to 
soft masses are negligible and the dominant source is radion mediation. The soft terms 
for bulk scalars in this case are (see, for example, (51IT7]) 



2R 



a-cf)7rkR 



sinh ((| 



cf)7rkR) 



(46) 



13 



and 
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- cf)2TrkR 
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-cf)2nkR _ 
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cf)2TTkR 



^)2nkR 



c!-)2TTkR 



■cf)2TTkR 




(47) 



It is well known that soft masses for brane-localized scalars cannot be induced by 
F"^ or F'^ at the tree level. Together with the Higgs mass degeneracy condition which 
holds in the NPT model, as discussed earlier, the Higgs mass parameters are then [6l|13] 



0, 



F'^kn 



\B„ 



(48) 



As explained in 13 , it is difficult to obtain realistic spectra with this condition. 



Finally, the gaugino mass is 



(49) 



As discussed in the previous Section, the gauge fields could be predominantly UV brane 
fields (if the theory is strongly coupled in the bulk, R is only moderately large, and there- 



fore the dominant contribution to the 4D gauge coupling in Eq. (41 ) comes from the UV 



brane term l/gfj^). In that case the radion- mediated gaugino mass will be suppressed, 
and since the brane contribution to is also small, no realistic phenomenology can 

be obtained. We are thus led to focus on the other case where R is large enough to 
overcome the NDA suppression of Eq. (42), RM^, ~ 0(100). Then the second term on 



the RHS of Eq. (41) contributes sizeably to l/g1; in the extreme case where the UV 
brane term can be neglected, we obtain 



M 



1/2 



F^ 
2R ' 



(50) 



comparable with the largest other radion-mediated soft masses. 



5 Parameterization and constraints 

As we showed in the previous section, if SUSY breaking is dominated by brane sources 
alone, then the gaugino masses will be relatively suppressed (to the extent that naive 
dimensional analysis is valid). In order to pass the constraints on chargino and gluino 
searches, the remaining soft terms would then have to be in the multi-TeV range. This 
scenario is clearly disfavoured on from the naturalness point of view, and would probably 
be impossible to probe at the LHC. Furthermore, large o-terms and soft masses as 
predicted NDA tend to lead to tachyonic sfermions in the low-energy spectrum. On these 
grounds we will dismiss the possibility that soft terms are induced by brane sources alone, 
and instead focus on the case where radion mediation gives a significant contribution. 

However, if the soft terms are exclusively generated by radion mediation and brane 
sources can be entirely neglected, we do not find realistic electroweak symmetry breaking. 
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In fact previous analyses 22 , 23 have shown that it is difficult to reconcile minimal 
radion-mediated scenarios with a GUT-scale degenerate Higgs mass matrix. While more 
refined scenarios [22u24] can give realistic TeV-scale physics, in our case the Higgs sector 



is subject to the even stronger condition Eq. (48), which turns out to be too restrictive. 



In short, radion mediation is necessary to provide sizeable gaugino masses, while 
brane sources are necessary to avoid vanishing Higgs soft masses. We will therefore 
study the general case where both sources of supersymmetry breaking are present. Their 
relative importance will evidently depend on the relative size of F^/2i? and F^/M,. It 
will also depend on tan/3, since we will determine the c-parameters from the Yukawa 
couplings which are fixed by tan /3 and the known fermion masses, and the soft masses 
depend on the Cj. 

The phenomenologically most problematic soft terms are the trilinear a-terms orig- 
inating from brane-localized SUSY breaking. Note that these are enhanced over the 



other soft terms according to the NDA estimate of Section 4.1 Large a-terms lead not 
only to large flavour violation but also to tachyonic third-generation sfermions; we find 
that to avoid these, the brane contribution should be subdominant with respect to the 
radion-mediated contribution, 

pZ pT 

Even with this condition satisfied, anarchic o-terms still tend to induce unacceptably 
large flavour violation — unless, again, the overall scale of SUSY breaking is unnaturally 



large. We therefore choose to set the brane trilinear term of Eq. (36) to zero, which is 
justified if the SUSY breaking fields on the brane are charged under some symmetry (for 
a comprehensive discussion see e.g. (71125]). 



To now investigate supersymmetric flavour violation, we decompose the brane- 
induced soft mass matrices into a hierarchical part which depends on the c-parameters 
and an 0{1) part: 



"^xjj, brane ~ 



F 



z 



2 



4A;^^ iX = U,D,Q,E,L), (52) 



where \"]^ are dimensionless hermitian matrices with 0(1) entries. Recall that k^^ was 
defined as 

~ M / \ 1/2 ' ^ ' 



with and Y^ as in Eqns. ^ and ([s]). Similarly, we write the Yukawa matrices yij as 

« 2 fU fQ 

u _ "'^ •'^ •'j \U 
^ij A/f / \ 1/2 *i ' 

„ 2 fD fQ 
D _ ^ ■'j \D 

~ M / A 1/2 ' (54) 



M / \l/2 

(^Yf Yf 

Here ^^j^'^ are dimensionless 0(1) matrices. Such a parameterization can be applied to 
any model which predicts the Yukawa matrices to be hierarchical, with power-suppressed 













1/2 


67r2 


fD fQ 






(YfY^) 


1/2 
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entries, up to a priori unknown anarchical C(l) coefficients. Examples include, besides 
our wave function localization scheme, also Froggatt-Nielsen type models. 

Even though the structure of the Yukawa matrices is dominated by the hierarchical 
part, it is important to take also the anarchical A-coefficients properly into account. An 
adequate framework for this is Bayesian statistics. We will treat the A-parameters as 
random variables with associated probability density functions (PDFs) /(A). The choice 
of PDF reflects our theoretical bias (for instance, that all matrix entries should be 0(1)), 
so this PDF constitutes a "prior" in the usual Bayesian vocabulary. The predictions of 
the model, such as masses or low-energy observables, will then also be PDFs. This 
approach will enable us to compute Bayesian credibility intervals, the equivalent of 
confidence intervals in the frequentist approach. 

We restrict ourselves to real couplings since the CP problem is not the subject of our 
study. It should however be kept in mind that, depending on the model, CP-violating 
observables may provide constraints which are just as stringent as those coming from 
FCNCs. We also neglect the neutrino sector, because it would induce additional model 
dependence. We focus on lepton flavour violation in the charged lepton sector, where 
the constraints are most stringent |26j. The wave- function suppression factors f^, and 
consequently the hierarchy structure, are determined by the six 5D bulk mass parameters 
cji and cjr. as described in Section [3] In addition to these, the following unknown 
0(1) flavour coefficients enter the analysis: Five symmetric 3x3 matrices X"^x for the 
squark and slepton masses in Eq. (|52|), a symmetric Yukawa coefficient matrix A^, and 



two unconstrained 3x3 Yukawa coefficient matrices A^'^ in Eq. (54); or a total of 54 
additional parameters. 

In the Standard Model subsector, we have the six c-parameters and the 24 Yukawa 
coefficients A]^'"'^. On the other hand, the experimental observables are the nine Stan- 
dard Model fermion masses and the three CKM angles. Simply setting all |Aj^'"'^| = 1 
and adjusting the six cji , cj^ , cji , cj-^ , cj-^ , cj-^ would not reproduce the Standard Model 
data with reasonable accuracy, so some deviation of the A"'*^'*^ parameters from unity is 
clearly needed. However, allowing all 24 to vary and attempting a full Bayesian 

analysis would require us to take into account the entire SUSY model, because of SUSY 
threshold corrections to the Yukawa couplings. This would be computationally very 
involved, and the results rather unwieldy. Instead we will take a simplified approach, 
which still allows us to extract the essential information. 

To quantify matrix anarchy, we allow the to vary independently within the 
range 

|Ai,-| G [1/^,^], (55) 

where ^ > 1 is a constant which is universal for all |Aij|. With logarithmic weighting, 
the prior PDF is 

/(log|Ai,|) = [/(-log^,log^), (56) 

U{a,b) being the uniform distribution on the interval [a,b]. We also allow the signs of 
the Xij to be independently ±1 (subject to certain restrictions; see below). As stated 
above, to satisfy the experimental constraints at some reasonable level of precision (e.g., 
1(T-3(t), a minimum ^ > 1 denoted by ^min is necessary. 

We also need to take experimental uncertainties into account. Each MSSM observ- 
able has an associated experimental PDF /ex, characterizing the uncertainty with which 
it is measured. On the other hand, any MSSM observable can be expressed in terms of 
our model parameters Xij, which for any given ^ defines its theoretical PDF f^. The 
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GUT-scale parameter 


tan /3 = 1.3 


tan /3 = 10 


tan/3 = 50 


yt 


4 (75%) 


0.48 (4.2%) 


0.51 (5.9%) 


yt 


0.0113(89%) 


0.051 (3.9%) 


0.37(5.4%) 


Vt 


0.0114(2.6%) 


0.07 (4.3%) 


0.51 (7.8%) 


yu/yc 


0.0027(22%) 


yd/ys 


0.051 (14%) 


ye/y^i 


0.0048 (4.2%) 


yc/yt 


0.0009(111%) 


0.0025 (8%) 


0.0023(8.7%) 


ys/yb 


0.014(29%) 


0.019(11%) 


0.016(12%) 


yfi/yr 


0.059 (3.4%) 


0.059(3.4%) 


0.05 (4%) 


A 


0.56(61%) 


0.77(2.6%) 


0.72 (2.8%) 


A 


0.227(0.4%) 


\p + ^vl 


0.397(16%) 



Table 2: Mean values and propagated relative uncertainties of the experimental PDFs /cx asso- 
ciated to the GUT-scale Yukawa matrix eigenvalues and CKM matrix Wolfenstein parameters. 



These propagated values are extracted from 27 . 



combination of these two PDFs gives the total PDF for any given observable]^ 

For each independent experimental constraint (or equivalently, each independent 
observable), there is now a a characteristic value of denoted by ^* , above which 
the total PDF is dominated by fj^. This happens when /ex and fj^ have roughly the 
same width. We take the /ex to be normal distributions of variance a^^, with GUT-scale 
propagated errors from |27| ; explicitly, the mean values and relative uncertainties are 
given in Table [2] The .if* associated to each constraint is estimated to be log .if* ~ ciex- 
While this is only a rough estimate, it will be sufficient for our purposes. 

The constraints with the largest ^* correspond to those that are the "hardest to 
fulfill" (i.e. to observables that are the hardest to fit). For a given constraint, once 
.if > .if*, its width (dominated by fj^) increases with .if, which makes the constraint 
easier to fulfill, until it effectively decouples from the fit. In a fit involving n constraints 
and p parameters, whose best-fit point has likelihood Lbf, the quantity — logLer will 
decrease as =5f increases. The fit becomes perfect, with — logLer = 0, once only p of 
the n constraints remain. This is illustrated in Figure [T] 

These observations can be used to determine =Sfmin as follows. We first perform a 
fit involving the six c-parameters, and the seven constraints with largest ^* . This gives 
a best-fit point with likelihood Lbf < 1- Using the PDF of the seventh constraint, we 
can then deduce the value of .if necessary to increase Lbf until the required level of 
precision is reached. This value of .if gives .ifmm, provided that the constraints which 
were previously not taken into account can also be satisfied with this .if. If one of these 
constraints is not satisfied, we repeat the procedure, including this additional constraint 
in the fit. 

The estimation of ^min depends on the precision required to fit the Standard Model 

^In the present situation, the PDFs are combined simply by convolution. For example, the top Yukawa 
coupling is given by yt = | A3^|e~^'^'°^''^T3 "'^Z^' to leading order (cf. Appendix [b| , hence —nkR{cj-^ — 
1/2) = (logyt - log |A|'3|)/2. The PDF of logy* - log [Xi^^l is the convolution of the logj/t and log jA^'al 
PDFs, in other words, of /ex(logi;t) and fse ~ [/(— log.jSf , log.5f ). The Xij can be regarded as nuisance 
parameters, the parameter of interest being ■ For the other couplings, the are more complicated 
because the dependence on the Xij is more involved (see Appendix IbI) . 
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- log 




Figure 1: Logarithm of the hkehhood of the best fit point, — logLer as a function of for a fit 
involving n constraints and p parameters. The constraints weaken and eventually decouple as ^ 
increases. The fit becomes perfect, i.e. — logLsF = 0, once only p of the n constraints remain. 

values (masses and CKM mixing angles) . It also depends on the running of the Yukawa 
couplings, which in turn depends on tan /3 and on threshold corrections. Finally, it 
depends on the efficiency of finding the best fit point. We find 

if„,i„ = 1.2 - 1.5. (57) 

The value ^min = 1-2 is obtained for a fit of 3a precision, while ^min = 1-5 corresponds 
to a more conservative la fit. It is reasonable to also impose an upper bound on our 
starting point being that the As should all be 0(1), but this bound is of course much 
less rigorously defined. In the analysis of Section |6j ^ will be allowed to vary within 
the range 

1.2 < ^ < 3. (58) 

Flavour matrix anarchy being an essential ingredient in our framework, ^ can in a 
sense be regarded as a model parameter which measures the allowed deviation from the 
superimposed hierarchical structure. 

Finally, to determine the c-parameters we use the constraints with largest * , which 
turn out to be the quark Yukawa couplings. The c-j-. are then deduced from the expres- 



sions given in Appendix pi subsection B.2, The cj^. , on the other hand, are not hier- 



archical, cjr^ ~ cj^2 — so they cannot be determined from the analytic expressions 
for the down- type Yukawa couplings. We therefore set cj-j = cj^^ = cj^ = cj^, with cjr 
determined by the bottom Yukawa coupling. We have checked that this choice does not 
sensitively infiuence the mass spectrum, the mixings, or the rates of flavour violating 
processes. 

We close this section with some remarks on the signs of the Ajj. In the limit ^ — )• 1, 
such that all Xij are ±1, some sign combinations give rise to accidental cancellations when 
rotating to the mass eigenstate basis, by which one or two of the fermion masses vanish 
exactly. For instance, if all Xij = +1 then there is only one mass eigenstate with nonzero 
mass. When ultimately taking ^ > 1, the previously vanishing masses will acquire 
widespread, ^-dependent PDFs. For greater predictivity we thus restrict our analysis 
to sign combinations which are non-singular in the limit ^ — )• 1. 

Furthermore, we expect the and Yukawa matrices to obey GUT relations 
at the level of signs. This is because all Yukawa couplings originate on the IR brane, 
where SU(6) is broken only spontaneously. Higher-dimensional operators involving some 
powers of (S) can lead to violation of the GUT relations, but the leading contributions 
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which determine the signs are given by the SU(6)-symmetric second term in Eq. (11), 
which imphes sign ( A^- j = sign ( Xj - j . The same argument holds for soft masses in the 



T and J- sector, that is sign 



m, 



'Qij 



Sign I 



"uij 



sign I m| 



and sign ( mlij 



sign I m-j 



"Dij 



6 Results 

We are finally in a position to describe our numerical analysis. The dominant flavour 
constraint comes from the lepton sector, more precisely from BR(;U — )• 67), as is the case 
for many other SUSY GUT models - see for instance ^28j. The leading contribution to 
/i — )• 67 contains tan /3-enhanced pieces 29 , so that in the mass-insertion approximation 



the branching ratio can roughly be estimated as 



BR(/i 67) 



a 1 







2' 













(59) 



where m-s is a typical sparticle mass, and m|, ml and are in the super- CKM basis g 
For our purposes, however, this estimate is not precise enough, and the branching ratio 
is instead determined by a full one-loop calculation 30 . The lepton sector comprises 21 



relevant A-parameters, 9 from the Yukawa matrix and 6 each from the soft mass matrices 
m| and . The soft terms in the quark sector enter only at higher loop order, so these 
are all the A-parameters which need to be fixed. 

Despite the fact that ^ should be larger than ^min in order to fit the Standard 
Model, we initially set .if = 1, so the only variables are the signs of the Xij = ±1. There 
are 3-2^^ physically inequivalent relevant sign combinations: By field redefinitions one 
can choose five of the signs in the Yukawa coefficients A^- to be positive, and of the 
remaining 16 combinations, only 6 lead to non- vanishing fermion masses for all three 

generations. Six independent signs in each X-^ and X-^ can be chosen independently, 
hence we have 6 • 2^ • 2^ = 3 • 2^3 ^ 24576 combinations. We first scan over these sign 
combinations while keeping ^ = 1, and subsequently allow for ^ > 1. 

For any given tan/3, the c-parameters are now determined as described in Section [5] 
For any given scale of radion- mediated SUSY breaking /2R and of brane-source 
mediated SUSY breaking /M^,, the soft terms at the GUT scale are calculated from 
Eqns. (45)-(48), (50), and (52). The sparticle mass spectrum, mixing matrices, and 
low-energy fiavour observables are computed using SPhenoS ^31j, appropriately modified 
to handle, in particular, the DHMM condition Eq. ( |44[ ). 

Figure [2] shows some sample distributions for BR(/x — )• 67) at ^ = 1 for various 
tan /?. At low tan /3, the dominant contribution to BR(/x — )• 67) comes from the trilinear 
term . By assumption, the trilinears are induced only by radion mediation and so 
do not depend on the A-parameters. This explains why in the left panel of Fig. [2] 
the value of BR(|U — )• 67) is hardly sensitive to the sign combination. By contrast, at 
large tan/3, and the brane soft masses w-lbranc ""^Lbrane (defined in Eq. (52)) give 



comparable contributions to BR(// — )• 67). Since these depend on the X-^ and A^^^ 
distributions span a much wider range at large tan /3. 



the 



■^In this formula the first two terms reflect the contributions of diagrams 1(c) and 2 of Reference 
while the last term corresponds to a single LR flavour-changing mass insertion. 
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Figure 2: Distributions of BR(/x — > 67) for all permissible sign combinations in the A matrices at 
if = 1. The model parameters are F'^ /2R = F^/M* = 1500 GeV, and tan/3 = 5, 10, 30 (from 
left to right). The current experimental bound is BR(/i ej) < 1.2 • 10^^"'^ 
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The arrows 

indicate the positions of two benchmark sign combinations SI and S2 used later in our analysis. 



To also allow for ^ > 1, and hence for a good Standard Model fit, we choose two 
benchmark sign combinations in the A matrices, which we denote by SI and S2 (indicated 
by arrows in Fig. [2]). Explicitly, we have 




o2 




A"^^ = + + + , A'"^ = + + + (60) 




for SI, and 






A"^'^ = + + + , X< = I + + + \ . (61) 



for 82. As is evident from Fig.[2j with the sign combination S2 the /i — >• 67 decay rate is 
suppressed even for large tan /3, whereas the SI point exhibits increasingly large flavour 
violation. We then scan several values of tan/3, /2R, and F^ /M^: for these two sign 
combinations, and allow the \Xij\ to deviate from 1. The values of the \Xij\ are drawn 



from their prior defined by Eq. (56) 



Figure [3] shows a typical example for the sign combination S2 and F'^/2R = F^/M^, 
such that (i.e. radion mediation) dominates at low tan/3, and brane soft masses 
dominate at large tan/3. In this case, BR(/x — )• 67) is suppressed at large tan/3. One 
can also see that the width of the PDF increases with ^ as expected. Note, moreover, 
that at large tan /3, where the brane contribution dominates, the mean value of the PDF 
is shifted towards larger values when increasing The reason for this ^-dependent 
shift is that certain cancellations occur between the brane contributions for this sign 
combination in the limit .jSf — )• 1. At low tan/3, where radion mediation dominates, this 
effect is less important. 

Regardless of the chosen sign combination, when radion mediation dominates, the 
BR(/i — )• 67) constraint is weakened by at least an order of magnitude, depending on .if. 
By contrast, when the brane-induced soft terms dominate, there are large contributions 
to BR(/i 67). 

In Figs. [4] and [5] we show the lower bounds on BR(/x — )• 67), given by 95% Bayesian 
credibility intervals (BCIs), for various .if ^ J^^mn- Figure |4j is for combination SI, while 
Fig.[5]is for combination S2. For any given value of the red regions to the right of the 
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Figure 3: Probability density functions of BR(/i — > ej) for the sign combination S2, with 
F'^/2R = F^/M., = 1500 GeV, and tan/3 = 5, 10, 30 from left to right. The red, green, 
blue lines correspond to .5f = 1.2, 1.5, 3 respectively, with the vertical black lines showing the 
values at ^ = 1. Radion mediation dominates at tan/3 = 5, while the brane source dominates 
at tan /3 = 30. The PDFs are normalized to have the same maximum. 



corresponding ^-line pass the experimental constraints]^ In the left panel, where radion 
mediation dominates, the bound on F-^/2i2 is weakened by a factor 2-3 when increasing 
. Since we have M1/2 ~ F'^ /2R, the gluino mass is typically rrig ~ 2 {F'^ /2R) due to 
RG running. Moreover, squark masses are dominated by the gluino contribution to their 
RGEs, so we have rriq ~ rrig. Therefore, the weakening of the BR(/i — t- 67) bound opens 
a part of the parameter space which is relevant for the production of SUSY particles at 
the LHC. 

We conclude from Figs. |4] and [5] that there are regions of the parameter space which 



pass the flavour constraints, and where SUSY particles are within discovery reach [33 34 
at the LHC. ATLAS and CMS searches with about 1 fb^^ of data at 7 TeV already put 
lower limits on gluino and squark masses of roughly rrig^q > 1 TeV for niq ~ r 
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36 . It is worthwhile noting that in our case the BR(;U — )■ 67) constraint forces the 
SUSY spectrum to be heavy, generically beyond the current LHC limits. In that sense 
the most severe constraints on our class of models still originate from flavour precision 
experiments, rather than from direct superpartner searches. 

At this point a comment is in order concerning the effects of subdominant flavour 
constraints. The next-to-dominant constraint turns out to be BR(/x — )• 3e), which is 
strongly correlated to BR(^ — t- 67), such that a weakening of the latter weakens also the 
former. We do not show this or other subdominant constraints here, since little would 
be gained by taking them into account. 

Production cross sections at the LHC are very similar to those of the mSUGRA case 
with rriq « rUg, and can be characterized by the gluino-squark mass scale, see e.g. [37||38j . 
At rrig^q « 1, 2 and 3 TeV, the overall SUSY cross section is of the order of 1 pb, 10 fb 
and 1 fb, respectively. For rrig^q ~ 1 TeV, the cross section is dominated by gluino- 
gluino, squark-squark and gluino-squark production. For heavier masses, rrig^q ~ 2- 
3 TeV, squark-squark and electroweak ino-ino (mainly xf X2) production dominate, 
while gluino production becomes negligible. However, interesting LHC signatures arise 
from the slepton mass patterns and mixings, leading to chargino/neutralino and slepton 
decays that are specific to the setup studied here. 

Let us therefore next discuss details of the spectrum, of mixings, and of LHC phe- 



^We consider nih > 114 GeV and nih > 111 GeV to account for the ~ 3 GeV theoretical uncertainty 
on the light Higgs mass. 
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Figure 4: The dominant constraints in the {F'^ /2R, F^ /M^,) plane for the sign combination 
SI, tan/3 = 5 (left) and tan/3 = 30 (right). F^/2i? and F^/AF are in GeV units. The 
black lines on the bottom left are nih = 111 GeV and nih = 114 GeV isolines. The blue lines 
show the BR(/i ej) constraint for ^ = 1, while green lines show 95% BCIs of the same 
constraint for ^ — 1.2, 1.5, 3. The red regions, towards large F^ /2R, satisfy rn/j > 111 GeV 
and BR(/i — > 57) < 1.2 • 10^"'^^ for at least one value of ^ . In the white regions, a too large 
F^ /Mif leads to tachyonic sleptons. 
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Figure 5: Same as Fig. [4] but for the sign combination S2. 



nomenology. Two cases will be distinguished. First, we will treat the case where the 
effects of F^ /M^ in the scalar soft masses are negligible compared to those of F'^ /2R. 
In this scenario, which we call radion dominated, the only effect of brane-localized SUSY 
breaking is to provide nonzero Higgs soft masses. Since our model has a GUT-scale de- 



generate Higgs mass matrix, Eq. ( 44 ) , it is convenient to fix the Higgs soft masses by this 
condition, and to set i^^/M* to zero for the purposes of calculating sfermion soft terms. 
We have checked that a small F^/M^,, of the order of ~ 0.02 {F'^ /2R), is sufficient to 
generate the necessary Higgs soft masses but has negligible effect in the sfermion sector. 

Second, we will discuss the situation where the scalar soft masses receive contribu- 
tions both from F^ /M^, and from F'^ /2R, which we call mixed brane-radion scenario. 
F^ /M^: should be bounded from above, because it can induce tachyonic sleptons through 
RG running if it is too large, and because it can enhance BR(/x — )• ej) depending on the 
signs of the A-parameters. This can also be seen in Figs. |4] and [5} 
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Figure 6: Spectrum of SUSY particles, norrualized to the gluino mass, when radion mediation 
dominates. The hghtest state is almost completely a right selectron cr. The left-handed slepton 
masses increase with tan/3, due to the y'^ enhancement (see text). For large tan/3 and large Jif, 
the RGE effects of flavour anarchy cause the slepton mass PDFs to grow wider, except the one 
for mg^. 

6.1 Radion dominated scenario 

In the case where scalar soft terms are dominated by radion mediation, the effects of 
flavour matrix anarchy appear only in Yukawa couplings, not in the soft terms. In 
other words, the soft terms do not directly depend on the A-parameters. They will 
be sensitive to flavour matrix anarchy, i.e to only through RGE effects. The mass 
ordering depends mainly on tan f3, and is charted in Fig.[6j Here we denote the six slepton 
mass eigenstates /i..,6 by their dominant components, for example li ~ cr, I2 ~ fj-R, and 
so forth. The li ~ cr turns out to be the lightest state of the spectrum, followed by 
the lightest neutralino Xi and the I2 ~ fj-R- The masses of left-handed sleptons and 
sneutrinos increase with tan/3. This is because, as tan/3 grows, the charged lepton 
Yukawa couplings increase, so the c-parameters shrink and the corresponding soft 
terms are enhanced. 

At small tan /3, the PDFs are strongly peaked at the values shown in Fig. [6j At 
large tan/3 and large .if (about 1-5 < .if < 3), the RGE effects of flavour anarchy widen 
the slepton mass PDFs, again due to the enhancement at large tan /3. This effect 
is hierarchical: For the ruf PDF it is larger than for the PDF, while mg^ shows 
almost no sensitivity. Typically, for tan/3 > 30 and ^ > 1.5, this uncertainty is large 
enough to flip the mass ordering of [Ir and Xii and of fi and X2- Other mass orderings 
are conserved. 

A stable charged slepton as the lightest SUSY particle (LSP) is obviously excluded by 
cosmology. In particular it is not a suitable dark matter candidate. Dark matter could 
instead be composed of gravitinos G, which are the LSP in many models of warped 
supersymmetry (such as the one in Appendix [A|), or of axinos a. In that case the cr 
decays as — )• -l-LSP (LSP = G or a). If this decay occurs at the epoch of Big Bang 



Nucleosynthesis (BBN) 39 , T{eR) > 1 sec, it can alter the yield of light elements. This 



poses important constraints in particular on the gravitino LSP case. In case of an axino 



LSP, the eR decay is much faster 40 , so that BBN constraints can be evaded easily. A 
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detailed discussion of this issue is beyond the scope of this paper; in the following we 
simply assume that eji is in fact the next-to-LSP and that its abundance and lifetime are 
small enough to evade cosmological constraints. Note, however, that even in the axino 
LSP case with T(eij) ^ Isec, the appears as a heavy stable charged particle (41 42 



in collider experiments. For definiteness, we will refer to the eji as the "lightest massive 
particle" (LMP) in the following. 

Let us now turn to LHC cascade decays. Gluinos and squarks, if produced, decay 
as 5 ^ qqR,L, Qr ^ qXi (~ 100%), and ^ q'xt (~ 65%) or qx^ (~ 30%). The x?,2 
and xf decay further, and the decay chains end with the cr LMP and an electron. If 
m^^ < rn^o, the Xi can decay both as x? — ^ ^^^^ and x? — ^ Ar/^^- The relative rate 
between the two decays is dictated by the ratio of mass splittings: 

BR(x?^e|e^) ^(^lo-^kV ^g^) 

^J■R J 

The (Ir decays dominantly through the three-body mode fiR — >• e/je/x via a virtual Xi- 
The decay width is typically of 0{keV), so there is no displaced vertex. (In principle the 
fiR may also decay through the LFV mode fiR cr + Z ii kinematically allowed, but 
this is suppressed by a very small coupling.) This contrasts with universal scalar mass 
scenarios, where the lightest slepton is typically the fi. Observing an electron or a muon 
associated to the LMP instead of a r at the end of the decay chains would therefore be 
a hint for non- universality in scalar lepton masses. 

As mentioned, the cr LMP is stable inside the detector and behaves like a heavy 
muon. This can be triggered on in the muon chambers of the ATLAS and CMS experi- 
ments. The muon chambers also allow excellent track reconstruction and time of flight 
measurements with an accuracy of around 1 ns, which should allow to reconstruct the 
mass of the LMP with good precision [IT]. Moreover, the rate of energy loss through 
ionization (dE/dx) may be used to identify the LMP and measure its properties |43|. 
Given this striking signature, Drell-Yan production of cr, even with a low cross section, 



may be exploited. See 41 and references therein for more details. 



A further, ambitious idea is to use a stopper detector to observe the LMP late de- 
cay 44 . With a sufficient number of events, flavour violating decays could be observed, 



and eventually used to gain some information on m|. Since m| is hierarchical, it in- 
duces a hierarchical mixing between the right-handed sleptons. LFV processes are thus 
suppressed by powers of the typical hierarchy factor e~'^^^ = e. In particular, one has 
roughly 

BR(eR ^ /i LSP) /BR{eR e LSP) ~ . 

Given that e should be 0{W~^) to reproduce the SM flavour hierarchy, if the branching 
ratio can be measured to 1% or better, this would provide a rough test for our scenario. 

The features discussed above are generic for our class of models. Other aspects of 
LHC phenomenology depend on the precise mass ordering, which in turn depend on 
tan/3. In the following discussion we will therefore distinguish the cases of small and 
large tan/3. For concreteness we will use three representative scenarios, one with low 
tan/3, and two different configurations with large tan/3 (one featurimg m^o < m^^, 
the other m^o > mjxn)- The spectra for our canonical choice of F"^ /2R = 1.5 TeV 
are given in Table [3j Since these points lie at the edge of the LHC discovery reach 
(total cross sections < 1 fb), we also provide in Table |4] an analogous set of points for 
F^/2fl = 1 TeV, which is more interesting for LHC studies. The complete SLHA files. 



including mass matrices and branching ratios, can be obtained from 45 . 
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For these benchmark points, we have also checked several observables from quark 
flavour violation using SUSY_FLAVOR. Due to the large theoretical uncertainty, the con- 
straints from flavour violation in the quark sector are all satisfied, so lepton flavour 
violation provides the strongest constraint as anticipated. 

In these points, the GUT-scale soft mass has been set to a universal value for 
simplicity, see end of Section [5} whereas generically the left-handed lepton soft masses 
may differ by 0{1) factors (they will, however, not be hierarchical since cj^^ ~ cj^^ ~ '^J^s)- 
A universal GUT-scale leads to small, RG-induced mass splittings between fiL and 
ei of typically 0(0.1%), and somewhat larger mass splittings between eL/fiL and fi 
of 0(5%). This is relevant because, at low tan/3, the X2 Xi are mostly wino and 
thus decay mostly into the left-handed sleptons ei, Ml, ^l- The branching ratios of LFV 
decays are at most ~ 10~^ in that case, and therefore irrelevant for LHC phenomenology. 

The main features can be summarized as follows: 



• The lightest Higgs mass myfi is always lifted by stop loop corrections to take a 
value closely above the LEP bound, as is typical for the MSSM with relatively 
heavy stops. 

• The xl decays as xl ^ '^f'^'^ (~ 20-30%) or xl ^ AiM'^/ef e=F (~ 10-15% each); 
the rest goes into i>jfj. The sleptons subsequently decay into Xi + s./ fi/r, followed 
by the Xi decay to the LMP, x? ^r^^- The resulting signature is X2 ~^ ^f^T^^Ry 
i.e. same flavour opposite sign (SFOS) dileptons, plus an electron, plus the LMP 
which behaves like a heavy muon. There is no Ej^^^^. 

• The xf cascade decays via a charged slepton II or sneutrino D into x^t^v^ (~ 50- 
60%) or into x\^^^e/ 20-25% each). The decay chain gives rise to dilepton 
signatures of e-|- (e/zx/r) with uncorrelated charges, plus the LMP, plus E^^^^ from 
the vs. This decay can be combined with the X2 decay on the other branch. 

• The masses of the sparticles appearing in the decay chains may be determined 
from kinematic distributions. The simplest observable is the endpoint of the SFOS 
dilepton invariant-mass distribution, MJ^^^ = m^o(l — m?/m|o)"'^''^(l — '"^|o/"i^)^'^^ 
from the X2 decayj^ For example, for point A', the Mn endpoints are M™^^ = 
308.28 GeV, M/^*^ = 308.61 GeV and M^'^^^ = 373.28 GeV. This may be used to 
obtain information on the masses of the three left-handed sleptons. 

At the LHC, it will most likely be not possible to reconstruct the GUT scale 
parameters through a bottom-up evolution. However, at the one-loop level, m^^ — 
mg^ is RG invariant, and the running of mf^ — mg^ depends on only a single 
parameter combination X-j- = 2\yr\'^{m'\j^ + mf^ + "^fj?) + 2|ot-P (cf. [47j). By 
combining the measurement of left-handed slepton masses with other information, 
such as the mass of coloured particles and the limits on (SUSY) LFV processes, 
one could at least carry out a hypothesis test on the structure of . 



At large tan/3, the situation is quite different. The left-handed sleptons and the 
sneutrinos are heavier, and the charged lepton Yukawa couplings are enhanced. The 
latter induces two effects through the RGEs. On the one hand, as explained above, the 
mass PDFs become much wider and either ordering, m^o > or m^o < m/ifl, can 

now occur. On the other hand, LFV processes in the — t sector may be sufficiently 

*Note however the ambiguity in the channel with electrons, X2 eje^ — >■ + 3e. 
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Figure 7: PDFs of x?-/ii? mass difference, for i^^/2i? = 1500 GeV and tan/3 = 30. The red, 
green and blue lines are for S£ = 1.2, 1.5 and 3, respectively. 



large to be observed. Since m^^ > "m^o, the sleptons relevant for LHC phenomenology 
are now the flu, cr and fi. The X2 decay predominantly into the /i" and, if 

kinematically allowed, into the fi, because it contains left-handed components induced 
by left-right mixing. The points B and C (B' and C) in Table |3] Q are representative 
examples. 

We will first discuss the leading decays, and then study LFV processes. Those features 
which are unrelated to LFV can be summarized as follows. 



Also in these benchmark points, the lightest Higgs mass m^jO is always around 120 
GeV, closely above the LEP bound. 

The Xi decays into e^e^ or, if kinematically allowed, also into fL^fj,"^ . If the muonic 
channel is open, its branching ratio strongly depends on the mass splittings, cf. 



Eq. (62). In this case the fiR decays predominantly through a three-body mode 
jln — eReiJ, as before. The signs of e and n are a priori not correlated, but they 
are related to the sign of cr, which may be measured, and the sign of the parent 
fljl. (If rufi^ > m^o, ji^ — )■ followed by Xi^^e^, giving the same final state.) 

The PDFs of m^o — m^^ are shown in Fig. jTj for different values of 

The X2 decays predominantly to a Higgs from X2 ~^ (~ 70%), if nif-^ > m^o. 
This is the case for points B, B' and C. If m^j < m^o (point C), the leading X2 

decays give a r+r" pair from X2 ~^ ^i^t^ ^ '^1^ ~^ 'T^Xi (~ 85%), while X2 ~^ Xi^" 
has ~ 10% BR. With the Xi decay, the whole chain produces an OS ditau or a 
plus the LMP, plus e 01 e + fj,^ + fi^ , depending on the fj-R/xi mass ordering and 
splitting. The final e ox e + fi^ + should be rather soft compared to the taus. 

Analogously, the decays mainly either through Xi ~^ Xi^^^ (~ 65-70%), or 
through Xi —7- f^^T — ^ Xi'T^'^T (~ 85%) if kinematically allowed. The latter chain 
gives one hard r plus EJp^^^, plus the LMP, plus e or e + + fi^ depending on 
the IIr/xi mass ordering and splitting. The subleading decay xf ~^ P'%'^ij. 
have 10-20% branching ratio and gives one hard fi plus E^^^^, plus e and the LMP 
if rrifij^ > m^o , or one hard fi plus E'™"^'^ , plus the LMP and e + fj, from the fiR 
three-body decay if m^^ < m^o . 

The invariant-mass distribution of the SFOS ditau can be used to determine the 
mass of the fi. For measuring the fiR mass, one needs to rely on the analysis of 
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Figure 8: PDFs of branching ratios relevant for LFV in the ji — t sector. The parameters are 
F'^ /2R = 1500 GeV and tan/3 = 30. The upper row shows ID PDFs for BR(n XiM)> 
BR(x2 — ^ A-R''') BR(x2 P-RlJ-)- The lower row shows 95% EC contours of the joint 2D 
PDFs of the these branching ratios with BR(t — > /i7). In all plots, the red, green and blue lines 
correspond to ^ = 1.2, 1.5 and 3, respectively. 



chargino decays, or on LFV processes, in which jlji appears as an intermediate 
decay product (see below). Knowledge of the masses, or mass splittings, of the 
fi, fiR and eji would now permit to obtain information on m|. As before, this 
information combined with other measurements would permit to carry out a hy- 
pothesis test, this time on the structure of m|. In particular, one can check wether 
the hierarchical factors e can be in agreement with values 0(10"^) necessary to 
reproduce the SM flavour hierarchy. 



Let us now turn to lepton flavour violation. As mentioned, LFV can be sizeable in 
the — T sector. The most interesting decays are those involving ^3 ~ fi and I2 — fiR- 
Figure [s] shows, in the top row, the PDFs for BR(fi — x^ji) and BR(x2 I^rt), and 
for comparison also for BR(x2 ft-RfJ-)- The bottom row shows the correlation of these 
branching ratios with BR(r — )■ ^^7). At the chosen scale, F'^/2R = 1500 GeV, the 
current experimental bound of BR(t — t- fi^y) < 6.8 x 10~^ is always satisfied. We observe 
that BR(r — t- fi'y) scales roughly as (F^/2i2)'' ' *■ Given the strong correlation, if 
LFV is observed in slepton or neutralino decays, this leads to a prediction for r — t- //7. 
Our main observations regarding LFV processes at the LHC are: 



• While the fi decays mainly to Xv^ can also have a LFV decay fi — )• Xi/^- 
The rate of the LFV decay peaks around 1% but can go up to 10% or more for 
.if = 3, see the top-left plot in Fig. [sj In the X2 decay chain, the ditau is then 
replaced by a /i^r^ pair, potentially giving rise to an interesting flavour structure 
in kinematic distributions. Kinematic edges with flavour splitting and mixing 
have very recently been studied in 148] (see also [49]). In the xf decay chain, the 
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logi„BR(x^^x!^'J) 

Figure 9: PDFs for cascade decays at F^/2i? = 1500 GeV, tan /3 = 30 and ^ = 3; dash-dotted 
line: BR(x§ ^ XiT=^t=f), dashed hne: BR(x^ ^ xV'^Ai'^): solid line: BR(x^ ->■ x?t±^=f). 

single T is replaced by a single ^, which must be separated from the non-LFV 
xf — )• p.Ri'fj, by kinematics. For example, one may exploit subsystem transverse- 
mass distributions [50] of, e.g., a lepton associated to the 'upstream' jet originating 
from ql cascade decays. This may permit to disentangle decay chains involving 
different slepton mass-eigenstates. 

• LEV can also occur directly in the X2 decays. Indeed the decay X2 ~^ i^R'^ can 
have a branching ratio of up to 0(10%) for .if = 3, larger than the X2 ~^ f^Rf^ rate, 
cf. the middle and top-right plots in Fig. [sj The decay chain then is X2 ~^ P-R^ 
XiMt) — >• eRtjiT with the x? being on- or off-shell depending on the P-r/xi mass 
ordering. 

The PDFs of the branching ratios of X2 ~^ Xi'^^t^^ X2 ~^ Xi^^f^^ ^-^id X2 ~^ 
Xi/i^A*^ are shown in Fig. [oj to illustrate the global rate of LFV expected in X2 
cascade decays. 



6.2 Mixed brane-radion scenario 

In the previous subsection we investigated the case /M^ ^ /2R, such that the 
scalar soft terms were dominated by radion mediation. We now consider the case where 
scalar soft terms receive non- negligible contributions from the brane source. The brane 
source contributions introduce a large uncertainty directly in the GUT-scale scalar soft 

terms. In addition to the \^-°'^ , the also become relevant. The phenomenology will 
eventually depend on the signs and magnitudes of all these A-parameters. Nevertheless 
we can still identify some generic features. 

Since the case of sizeable F^ /M^, is much more constrained by BR(^ — )■ 67) at large 
tan/3 (see Figs. [2| [4] and [S]), we will focus on the small tan/3 scenario. In this case, the 
left-right mixing is negligible, because the a-terms are not large. It turns out that the 
right-right mixing is hierarchical, as in the radion-dominated case. We therefore call the 
right-handed sleptons cr, /xr, and tr. On the other hand, the left-left mixing can now 
be very large, so we denote the left-handed sleptons by Ili, Il2, and II^. 

The masses of the right-handed sleptons can now span a wide range. All possible 
mass orderings with respect to the lightest neutralino can appear: < m^o < m^^, 

me^,/i^ < Tn^o, or ni^o < rrie^^ftn- This last possibility is particularly interesting since. 
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Figure 10: Probability of having a Xi LSP as a function of {F^ /M^,)/ {F'^ /2R). The red, green, 
blue lines correspond to ^ = 1.2, 1.5, 3, respectively. The other parameters are tan/3 = 5, 
F^/2R = 1.5 TeV. 



unlike the radion-dominated case, it features a Xi LSP which is a viable dark matter 
candidate. Matrix anarchy again plays a crucial role in realizing this possibility: With 
^ = 1, the RG invariant 

^ = '^Hu - "T-L + Tr(m^ - ml -2ml + ml + (63) 
would vanish due to exact SU(5) relations, and the LMP would then be the cr as in 

the radion-dominated case. However, with ^ > 1, the A^^-^ can be different from one 
another and induce a non-zero S. If S is sufficiently large and negative, the lightest 
slepton mass can be lifted above the neutralino mass. The probability of finding a 
neutralino LSP thus depends on as well as on the ratio {F^ /M^) / {F'^ /2R). This is 
illustrated in Figure 10, which shows the probability of finding a Xi LSP as a function 

2 

of ^ and {F'^ /M^)/{F'^ /2R) for a favourable sign combination of A-^^. This plot is for 
tan (3 = 5 and F^^ /2R = 1500 GeV, but the result is fairly insensitive to the SUSY scale. 

We proceed to discuss the left-handed slepton masses. At small tan /3, the soft masses 
ml are suppressed with respect to the gaugino masses (as explained in subsection 6.1). 
The RG running of the masses of the left-handed sleptons is therefore dominated by the 
gaugino masses and the S parameter. For {F^ /M^)/{F'^ /2R) < 1.5, we obtain the mass 
ordering 

m^o < mj^ < m^o . (64) 

This property is particularly interesting for LHC phenomenologyj^ 

For the remaining discussion, we fix tan /3 = 5, F^ /2R = 1500 GeV (i.e. mg ~ 
3 TeV). It turns out that, even for small F^ /M^^, the LFV effects in the SUSY decays 
can be large, particularly in the e— // sector, while still satisfying the current BR(;U — )• 67) 
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bound. The details depend on the signs in both \f- and ■ We will focus on the 
X2 ~^ Xihlj decays. Following 51 , we define the observable Kij as 

K- . = BR(x^ j^q^^ 

BR(x§^/,H^X?) + BR(x^^/f/Jx?)' ^ ^ 



^ For {F^ /M,)/{F'^ /2R) > 1.5 with ^ = 3, the S parameter of Eq. (|63| can be sufficiently large to 



make the left-handed sleptons lighter than the right-handed ones. This is why, in Figure 10 P{xi LSP) 
decreases above a value of {F^ / M*) / {F'^ /2R) > 1.5. Due to D-term splitting between sneutrino and 
charged slepton masses, a sneutrino becomes the lightest SUSY particle of the spectrum. Left-handed 
sneutrino dark matter is strongly constrained by direct detection and cosmology, so one would again have 
to assume that the actual LSP is the gravitino or an axino. We will not pursue this case any further. 
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Figure 11: Top left: PDF of BR(x2 ^ hhxi)- Plain, dashed and dash-dotted hnes correspond 
to Uli = ee, /i/^, TT respectively. Top center: PDF of BR(x2 ~^ hljXi)- Plain, dashed and 
dash-dotted lines correspond to Ulj — e/i, er, /ir respectively. Top right: Plain, dashed and 
dash-dotted lines correspond to i^e/^j K^t, K^^ respectively. Bottom: 95% BC region of the 
joint PDF of BR(/i ^ 67) with BR(x^ ^ e^Xi) (left) and K^^, (right). 

All of these distributions are for a fixed sign combination at F^/2i? = 1500 GeV, F^/Af* = 
50 GeV and tan /3 = 5. The chosen sign combination represents the typical behaviour of SUSY 
BRs, while it is somewhat favourable regarding the BR(/x — >• 67) constraint. 



to quantify the rate of LFV in the X2 ~^ ^f^JXi decays. Figures |lT| and [l2| show the PDFs 
of X2 Xihlj decay branching ratios, with F^/M^ = 50 GeV and FyM^ = 1500 GeV, 
respectively. For /M^, = 50 GeV, it is in particular LFV in the e — fi sector that can 
be large enough to give observable effects at the LHC, while at F^/M^ = 1500 GeV, 
LFV in all three sectors can be sizeable (though LFV still tends to be largest in the 



e — /U sector). The distributions shown in Fig. 11 and 12 represent the typical behaviour 



as far as BRs of SUSY particles are concerned, but the sign combination is chosen such 
that it slightly favours small BR(^ — )• 67). Note that, as opposed to Figure [sj there is 
no strong correlation between BR(x2 ~^ c^/i^Xi) ^-^id BR(^ — )• 67). 



Let us now discuss possible LHC signatures. Collider signatures will again depend 
on the various possible mass orderings of the Xi with respect to the and /i/j. There 
are the following possibilities: 



If m^o < rriej^^fij^, then the x? is stable and a dark matter candidate (barring the G 
and a LSP options). In the detector, the decay of X2 produces OS dileptons plus 

ipmiss 

If mg^ < m^o < m^^ , or if m^o 
eeR. The X2 therefore decays as X2 ■ -i -j 
plus an electron, and the track of eji, without 



-rrifiji is sufficiently small, the Xi decays mainly to 
^ iflJecR. The signature will be OS dileptons 

miss 
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If me^,/i^ < m^o and |m/i^ — mg^jl is not too small, the Xi has a sizeable branching 
fraction to eeji and ^/i_R, and the next-to-lightest slepton decays to the lightest 
slepton within the detector. We will assume in the following that the lightest 
slepton is the cr. The fiji decays through the three-body decay fiji — )• efieji. The 
complete X2 decay chains are thus X2 ~^ ^t^JXi ~^ iflJ^^R or X2 ~^ ^t^JXi ~^ 
iflJlli^fi'^jeeR- The signature of the X2 decay will be OS dileptons plus e, and 
the track of e/j, without E^^^^ and possibly with additional OS dimuons. 

If the fin is guaranteed to decay within the detector, there is no ambiguity in the 
observation of LFV. One can detect LFV by simply counting leptons: These events 
will have an odd number of leptons (3 or 5). Flavour conserving events will have 
an odd number of e and an even number of fx and r, while LFV processes will 
have either an even number of e and an odd number of fi or r, or an odd number 
of e and an odd number of /U and r. The observable Kij, defined in Eq. (65) to 
quantify LFV, can be directly measured by lepton counting : 

^ ^ BR(x^ ^ lflJ^i[t^^l^^]eeR) 

Finally, if Im^^j — mg^l is sufficiently small, both sleptons are stable within the 
detector. The X2 decays are therefore either X2 ~^ ^t^JXi ~^ iflJ^^R or X2 ~^ 
ifljxi — )• lflJfJ,fLR. The signature of the X2 decay will be OS dileptons plus e or 
H, and the track of cr or /i/j, without E^^^^. LFV in the e — t or fi — t sectors 
may be observed by detecting a single r in these decays. The observables Ker and 
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(67) 



(68) 



Kfj^r can be inferred without ambiguity using 

^ ^ BR(x^ ^ T^e^ecR) 

BR(x^ ^ r±rTeeR) + BR(x^ ^ e±eTeeR) 
iV(r±e^e) 
~ Ar(r±rTe) + iV(e±eTe) ' 

^ ^ BR(x^ ^ T^ti^HfiR) 

BR(x^ ^ r±rTMAii) + BR(xO ^ m^/^^/^Ar) 

.¥(r^^^/i) 

On the other hand, in the e — /i sector, one cannot detect LFV by simple lepton 
counting since ee// or e/in combinations can be produced both in flavour conserving 
and in flavour violating channels. However, the flavour-conserving channels give 
same- flavour dileptons with opposite signs (SFOS), while in the flavour-violating 
channels all sign combinations of the SF dileptons appear with equal probability. 
To disentangle the two contributions, one should look out for same-sign dileptons, 
which can only appear through the flavour violating channels. Neglecting the 
flavour violating effects of the right-handed slepton sector, which are O{10~^) 
(more precisely 0{e^)) because of the hierarchical mixing, Ken given by 



2 X BR(x^ ^ e±e±/xej?) 



BR(x^ ^ e^e^een) + BR(x^ ^ ti^fi^een) - BR(x^ ^ H^H^esR) 

2 X N{e^e^fi) 



(69) 



iV(e±e=Fe) + N{iJ,^fj,^e) - N{iJ,^fj,^e) 

or 



2 X BR(x^ ^ f^^fi^efiR) 
i + I 

2 X Nifi^^^t 



BR(x^ ^ li^fi^fijlR) + BR(x^ ^ e^eTf^fiR) - BR(x^ ^ e±e±M/iii) 



(70) 



The additional term in the denominator supresses the contributions to BR(x2 ~^ 
e^e^ lifxR) [ BR(x2 I^^P'^eeR)] coming from flavour violation. This ensures that 
the denominator contains only flavour conserving contributions. 
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7 Conclusions 



We have studied flavour violation in supersymmetric models with a GUT-scale warped 
extra dimension. With matter fields located in the bulk, and the Higgs fields as well 
as the SUSY breaking hidden sector localized on the infrared brane, exponential wave 
function profiles can at the same time generate hierarchical fermion masses and mixings 
and somewhat suppress flavour changing neutral currents. 

However, we find that the constraints on FCNCs (in particular those on lepton 
flavour violation) are stringent enough to still rule out most generic models. For the 
concrete example of the holographic GUT model of NPT, several additional assump- 
tions on the hidden sector are necessary in order to obtain a realistic spectrum and 
evade the experimental bounds. More specifically, there should be contributions to the 
SUSY breaking soft terms both from the radion superfield and from the brane-localized 
hidden sector fields; and the brane- induced trilinear soft terms should be small or zero 
(which could be enforced by symmetry). With these assumptions, substantial regions 
of the parameter space can give rise to realistic sparticle mass spectra while avoiding 
unacceptably large lepton flavour violation. 

In these surviving regions of parameter space, the LHC phenomenology depends 
on whether the soft terms are predominantly induced by the radion, or whether the 
contributions from the radion and from the brane-localized hidden sector fields are com- 
parable. We have given an account of the expected mass spectra and LHC signatures in 
both cases. 

Generically, the bounds on lepton flavour violation, in particular BR(e — >• ^^7), force 
the spectrum to be heavy, with squark and gluino masses well above 1 TeV. One the 
one hand this might explain why no signal of SUSY has yet been observed at the LHC 
operating at ^/s = 7 TeV. On the other hand it means that, should the setup studied 
here be realized in Nature, it will require the high-luminosity run at 14 TeV to explore 
it. Moreover, a detailed experimental study of this scenario will most likely require 
precision measurements at even higher energy and/or luminosity (LHC upgrade) [53] , 
or at a multi-TeV e^e~ collider |54|. 
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Table 3: Sample spectra for F"^ /2R =1.5 TeV. Points A-C are representative for the radion- 
dominated scenario with low and high tan/9: point A has tan^ = 5, while points B and C have 
tan/3 = 30 with respectively m^o < m/i^ and m^o > mjXf^. Point D is an example of a mixed 
brane-radion scenario with a neutralino LSP. 
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Table 4: Sample spectra for /2R = 1 TeV. Points A'-C are representative for the radion- 
dominated scenario with low and high tan/3: point A' has tan/3 = 5, while points B' and C 
have tan/^ — 30 with respectively m^o < m^^j and m^o > m^^. Point D' is an example of a 
mixed brane-radion scenario with a neutralino LSP. All masses are in GeV units. Production 
cross sections (in fb) at the LHC with -^s — 14 TeV, computed with MadGraph 
given. 
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Appendix 



A Radius stabilization and SUSY breaking 



In this Appendix we sketch a warped 5D model (following |55[[56] ) in which the radion is 
stabilized, and in which both the radion superfield and some additional IR brane fields 
have non-vanishing F-term expectation values. It therefore provides a dynamical origin 
for the SUSY breaking field background on which we based our analysis. 

Consider first a warped extra dimension with two separate sectors: A super- Yang- 
Mills theory in the bulk and a super- Yang-Mills theory on the IR brane. In the infrared 
the degrees of freedom are a radion superfield T and two non-abelian gauge superfields 
which will undergo gaugino condensation. The strong-coupling scales of the SYM theo- 
ries are taken parametrically smaller than the KK scale. Later we will add an "uplifting" 
sector, consisting of a dynamical SUSY-breaking sector on the IR brane. 

In units of the 4D reduced Planck mass M4 = 2.4-10^'^ GeV, the effective La grangian 
after gaugino condensation can be written as 



C 



-3e 



+ / d 



Here the radion Kahler potential is 

K = -?. log 



Ml 
k 



1 



+ c] + h.c. 



(71) 



(72) 



and a, b, c are constants. While a and b come from bulk gaugino condensation, and are of 
order unity (or somewhat large since the theory is weakly coupled at the compactification 
scale), c comes from the IR brane gaugino condensate and is exponentially small. Note 
that our conventions differ from those of [55j by the sign of k and by a factor vr in the 
definition of the radion field. 



Eqns. ( 71 ) and ( 72 ) are written in a frame where the warp factor is unity in the IR 



and exponentially large on the UV brane. For consistency with our conventions in the 
main text, we perform a Weyl rescaling, which is a redefinition of the chiral compensator: 



TfcT, 



This gives 



C = J d'^e^ip (^-3e-^/3) + j d^Oif^ {ae 



with 



K 



-3 log 



Ml 



1 



,-kTT{T+T) 



(73) 

(74) 
(75) 



It is convenient to define the warp factor superfield uj by the holomorphic field redefinition 

uj = ^e-^^^ (76) 



(note that this does not just amount to undoing the Weyl rescaling of Eq. ( 73 ); the chiral 
compensator is always normalized such that (93) = 1 + while here we are choosing 

a diff'erent way of parameterizing the radion). The Lagrangian becomes 



C 



3M| 



d'^d (^99 -ujoj)+ I d^e {auj^+^if-" + cw^) + h.c. , 



(77) 
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where v = b/kir. This yields the F-terms 

F'^ = ^auu'+\ 

3M| 

and the scalar potential 

At large warp factors, i.e. small |a;|, the \F'^\'^ term in V is subdominant, and the 
potential is minimized at a finite value of u, 

3c 



U = — . (80) 

(3 + 1^)0 ^ ^ 

This ansatz is self-consistent because c is exponentially small. There is also a decom- 
pactification solution at w — )• 0, which is however of no interest for us. 

For the gravitino mass we find 

™i/.-w=»^(^)^ci^M«~iF-r (81) 

Returning to the old variables, the radion F-term is 

' - — ] . (82) 



2R 27rkR \ uj , 

juj vanishes to leading order, but the subleading terms turn out to be finite and are 
parametrically of the order O 

So far the vacuum energy density is negative, and the vacuum is a non- 
supersymmetric AdS minimum. This can be remedied by adding an additional sector 
which breaks supersymmetry dynamically on its own, in the rigid limit, thus providing 
a positive contribution to the cosmological constant. The details of such an "F-term 
uplift" have been worked out mainly in the context of effective field theories from type 
JIB flux compactifications (see e.g. (57^). Our main interest here is the relative impor- 
tance of the contributions to soft terms from the uplifting sector compared to the radion 
contributions. In our normalization, including SUSY breaking IR brane fields Z/ in this 
background as 

A£brane = ^ d^B e-'^''^^^\Zif ^ J d^O e'^^^^^W [Z i) (83) 
will give an additional contribution to the scalar potential, 

^g-2^fc/? ^|^Z,|2 ^ (84) 

To fine-tune the cosmological constant to zero, we thus need 

Therefore, in this particular model, the brane-localized contributions to the soft terms 
dominate over the contributions from the gravitational sector, i.e. the radion and com- 
pensator contributions. Since the warp factor is only moderately large in the scenarios 
we are considering in the main text, this model could still serve as an example for mixed 
brane-radion mediation. The gravitino is naturally the LSP, as is common in models of 
warped supersymmetry. 
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B Matching parameters to fermion masses and mixings 



B.l Preliminaries 

The superpotential of the MSSM contains the quark Yukawa terms 

W = y[j HuQW + yf^ HdQ'D^ . (86) 
Both Yukawa matrices and y° may be diagonahzed by bi- unitary transformations, 

y^^,^=UuLy"ulj^, yi,,^=UDLy^U^DR- (87) 

Here Uul is a unitary matrix which diagonahzes the Hermitian matrix y^y"\ 

UUL {y"y"^) 4l = ydiag^diag^ , (88) 

and IAdl is a unitary matrix which diagonahzes the Hermitian matrix y°y°\ 

UdL [y^y^^) 4l = ydiagydiag^ , (89) 
The CKM matrix Vckm of quark mixings is given by 

VcKM = UulU^dl . (90) 

and the physical quark masses are given by the matrix entries of y^j^g and y^iag multiphed 
by the appropriate Higgs expectation value. Three mixing angles, one phase, and six 
mass eigenvalues constitute the physical observables in the quark sector. 

If there is only one small suppression parameter e, as in the simple example in 
Section [3j the Yukawa matrices are 

y" = ■ (91) 

\ Afie Afge Afge / 

Here e ~ 0.1, and the A^'*^ of order unity. Note that to leading order in e, the structure 
of y"y"t and y'^y'^^ is similar, up to an overall factor: 

yUyU\ ^ ^11^3 ^^,2 ^ ydyd^ ^ ^2 I ^3 ^^,2 I , (92) 











y'= 1 




^22^ 


-^23^ 






'^32^ 


^33 



4 

rue 


Sue^ 


tue' 








t* 




Wu 



4 


Sue^ 


toe' 




Uoe' 


Voe 


tie' 


v^e 





Here we have defined 



Tu 
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+ ^11-^31 




Vu 


\U \U* , \U \ u* 2 
— ^23^33 "T ^22^^32 ^ 


+ ^21'^31 





and 



— \ u \ u* I \u \ u*2 I \u \ u*A 
u — A]^3A23 + e + A11A21 e , 

, «. = |A2^3|2 + |A2^2pe2 + |A2"i|V, (93) 

, ^«. = |A3^3P + |A3"2|V + |A3-1|V, 



Td — |Aii| -t- |Ai2l "T |Al3| , So — AX1A21 -t- AX2A22 i- AX3A23 , 

= A11A3* + Af2Af2 + Af3A3|, "Ud = |A2l|^ + |A22|^ + IA23P, (94) 

_ \D \D* I \D \ D* I \D \D* _ \ \D \2 j_\\D \2 j_\\D \2 

— A21A31 + A22A32 + A23A33 , WD — \A'jl\ + IA32I + IA33I . 
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B.2 Fermion masses 



A hermitian matrix of the form 



M 




has the eigenvalues 



(95) 



1^2 



W + 



uw 



w 



w 



-{uw-\vf) + \tne^ + Oie''), 



w 



\V\' 2 



\sw — tv 



*|2 



u 



H3 = - {{ru- |sp) - 



uw 
\sv — tu\ 



\v\ 



iuw-\v\^) e* + C>(e''), 



(96) 



uw 



\v\ 



Note that in the case of M = y^y'^\ all of the "minors" {ru — |sp), {sw — tv*), {sv — tu) 
and {uw — \v\^) are 0{e^) (unless there is some fine-tuned cancellation between the 
A|i). Their ratios arc therefore 0{1), and expressions such as {sw — tv*Y / {uw 



If P) or 



{sv — tu) /{uw — \v\ ) are 0{e ). In particular, the smallest eigenvalue /X3 = y„ is only 
generated at higher order, namely at 0{e^). 

For the case of real A, we thus find the Yukawa couplings 



yt = + 



(A2^3)' + (Af2)' 



22 



2 IA33I 

^32^231 ,2 



e' + 0{e% 



e' + 0{e'), 



Vu 

Vb 



Vd 



0{e% 
/^e + C»(e^), 



(97) 




««-^e^ + 0(e"), 

Wd 



In many grand-unified models, including our benchmark model, the matrix A|J- is sym- 
metric because it originates from a lOjlOjS/^ coupling (in SU(5) notation), so A23 = A32. 

For the more general case that there are three distinct c-paramctcrs in the up-type 
sector, one should parameterize the wave function suppression in a more general way, 
by allowing for several distinct suppression factors. Defining 



£. = g-7rfeii(|cri|-l/2) 



(98) 



the up-type Yukawa matrices have the structure 



/ 


(if 


ei £2 


£1 


\ 






(£3)" 


£3 










52 






(ear 


£3 




\ 


ei 
£3 


12 
£3 


1 


/ 



(99) 
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where we did not explicitly write any A factors. Assuming that es ^ e2 ^ ei, the 
Yukawa couplings of the first two generations are, to leading order, 



„ — \ \U \(^ \2 j_ (-^23)^ + (-^32)^ / n2 , 

- 1-^3311^3; H r, I X [/ 1 i^2j +•••, 



Vc 



2 l-^ssl 

I'^33'^22 ^^32^231 n2 , 

— — + 



(100) 



In the down- type sector, with the assumption cj^i c_f2 ~ cj-3, all that needs to be 
done is to set 



e = e 



-7rfcR(|c^,|-l/2) 



(101) 



in the last three of Eqns. (97). 



B.3 CKM matrix 



Again for the case of a single small parameter e, the matrix M of Eq. (95 ) is diagonalized 
by 

Ul=\ Ce l-i(|7p + |CP)e2 -je \+0{e'), (102) 
V 'i^e' Ye l-\\l?e^l 

where 

V ^ sw — tv* 
7=-, C = 2- 103 

Applying this to calculate Uul and Udl-, we find for the CKM matrix 
VcKM = ^ulI^dl = 

( l-U\Cn\' + IC.P - 2C.a) (C. - Cu)e - + f^^^^fp^) 

+ 0(e'). 



(104) 



Note that the leading terms in 'ju and C have a rather simple form. 



23 , /n/'^2\ /- ^12^33 ^13^32 , fn^.^2\ (-ir\K\ 
"^33 ^22^33 ^^23^^32 

while 7b and are fairly complicated when expressed in terms of the Xfj. 

In terms of our parameters, for all A real, the CKM mixing angles O12 and 623 are 
then approximately given by 

sin 6*12 = (Cd - Cu)e, sin 6*23 = (Id - 7u)e ■ (106) 
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